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ABSTRACT 

The Sunyaev-ZePdovich Effect (SZE) has been observed toward six massive galaxy clusters, at red- 
shifts 0.091<z<0.322 in the 86-102 GHz band with the Y. T. Lee Array for Microwave Background 
Anisotropy (AMiBA). We modify an iterative method, based on the isothermal /3-models, to derive 
the electron temperature T c , total mass M t , gas mass M gl and integrated Compton Y within r 2 5oo, 
from the AMiBA SZE data. Non-isothermal universal temperature profile (UTP) f3 models are also 
considered in this paper. These results are in good agreement with those deduced from other ob- 
servations. We also investigate the embedded scaling relations, due to the assumptions that have 
been made in the method we adopted, between these purely SZE-deduced T c , M t , M g and Y. Our 
results suggest that cluster properties may be measurable with SZE observations alone. However, the 
assumptions built into the pure-SZE method bias the results of scaling relation estimations and need 
further study. 

Subject headings: cosmology: observation — galaxies: clusters: — sunyaev-zeldovich effect: 



1. INTRODUCTION 

The Sunyaev-Zel'dovich Effect (SZE) is an useful tool 
for studies of galaxy clusters. This distortion of the Cos- 
mic Microwave Background (CMB) is caused by the in- 
verse Compton scattering by high energy electrons as 
the CM B propagates through the ho t plasma of galaxy 
clusters (jSunvaev fc Zerdovicld[T972T ). The SZE signal is 
essentially redshift independent, making it particularly 
useful for determining the evolution of large-scale struc- 
ture. 

For upcoming SZE cluster surveys (IRuhl et al.ll200i 
iFowlerl 120041: IKanekdl20ul IHo et all 12009ft . it is impor- 
tant to investigate the relations between SZE flux density 
and other cluster properties, such as mass, temperature, 
and gas fraction. By assuming that the evolution of clus- 
ters is dominated by self-similar gravitational processes, 
we can predict simple power law relations betwe en inte- 
grate d Compton Y and other cluster properties (|Kaiserl 
1986). Strong correlations between integrated SZE flux 
and the mass of clusters are also suggested by numer- 
ical simulat ions (|da Silva et all 12004: iMotl et all 120051 : 
iNagail [20061) . These relations imply the possibility of de- 
termining the masses and temperatures of clusters, and 
investigating cluster evolution at high redshift, with SZE 
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o bservation data a lone. 

Uov et all (|2001D and iBonamente et all (|2008fl demon- 
strated an iterative approach based on the isothermal (3 
model to estimate the values of electron temperature T e , 
total mass M t , gas mass M g , and Compton- Y from SZE 
data alone. In this paper, we seek to derive the same 
cluster properties from the AMiBA SZE measurements 
of six clusters. Due to the limited u — v space sam- 
pling, the AMiBA data do not provide useful constraints 
on the structural parameters, /? and r c , in a full iter- 
ative model fitting. Instead, we adopt (3 and r c from 
published X-ray fits and use a Markov Chain Monte- 
Carlo (MCMC) method to determine the cluster proper- 
ties (T e , M t , M g and Y). We also estimate these cluster- 
properties from AMiBA data with structural constraints 
from X-ray data using the non-isothermal u niversal tem- 
perature profile model (Hall man et al.ll2007l ). All quanti- 
ties arc integrated to spherical radius r25oo within which 
the mean over-density of the cluster is 2500 times the 
critical density at the cluster's redshift. We then investi- 
gate the scaling relations between these cluster properties 
derived from the SZE data, and identify correlations be- 
tween those properties t hat are induc e d by the iterative 
method. We note that iHuang et all (|2009) investigate 
the scaling relations between the values of Compton Y 
from AMiBA SZE data and other cluster properties from 
X-ray and other data. All results are in good agree- 
ment. However, we are concerned that there are embed- 
ded relations between the properties we derived using 
this method. Therefore, we also investigate the embed- 
ded scaling relations between SZE-dcrivcd properties as 
well 

We assume the large-scale structure of the Universe 
to be described by a flat ACDM model with fl m = 0.26, 
Q A = 0.74, and Hubble constant H = 72 kms" 1 Mpc" 1 , 
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correspon ding to the values obt ained using the WMAP 5- 
year data (|Dunklev et al.ll2009l ). All uncertainties quoted 
are at the 68% confidence level. 

2. DETERMINATION OF CLUSTER PROPERTIES 
2.1. AMiBA Observation of SZE 
AMiBA is a c oplanar interferometer (|Ho et al.l 120091 : 



iChen et all 12009ft . During 2007, it was operated with 7 
close-packed antennas of 60 cm in diameter, giving 21 
vector bas elines in u-v sp ace and a synthesized resolu- 
tion of 6' (|Ho et all 12009ft. The antenna s are mounted 
on a six- meter platform ( Koch et al.ll2009i ). which we ro- 
tate during the observations to provide better u-v cover- 
age. The observations of SZE clusters, the details about 
the transform of the data into calibrated visibilities, and 



the es timated cluster profiles are presented in Wu et al 



2009). Fu rther system ch e cks ar e discussed in lLin et al 



2009) and Nishioka lit~"aTl J2009). For other scientific re- 



sults deduced from AMiBA 2007 observatio n s, please re- 
fer to lHuang et al.l (12009ft: [Liu et al.l (12010ft : iKoch et al.1 
(|2010ft : lMomar et al.l (|201oh : lUmetsu et all (|2009ft 

2.2. Isothermal (3 modeling 

Because the u-v coverage is incomplete for a single SZE 
experiment, we can measure neither the accurate profile 
of a cluster nor its central surface brightness. There- 
fore we have chosen to assume an SZE cluster model 
and thus a surface brightness profile, so that a corre- 
sponding template in the u-v space can be fitted to the 
observed visibilities in order to estimate the underlying 
model p arameters. We consider a spherical iso thermal 
/3- model (|Cavaliere fc Fusco-Femianolll976. 1978ft . which 
expresses the electron number density profile as 



nJr) = n n 1 + — 



-3/3/2 



(1) 



where n C Q is the central electron number density, r is the 
radius from the cluster center, r c is the core radius, and 
(3 is the power-law index. 

Traditionally the SZE is characterized by the Compton 
y parameter, which is defined as the integration along the 
line of sight with given direction, 

y(n) = [ a T n e ^^dl. (2) 

Jo m cC 2 

Compton y is related to A/sze as 



A/sze = IcMByf{x,T c ) 



xe 



1 



(3) 



where x=hv /^bTcmb, ^cmb is the present 

CMB specific intensity, and f (x, T e ) = 

\x cot h(.T/2) - 4] [1 + 5 Te \(x, T c )\ (e.g., (LaRoaue et al.l 
12006ft. J re i(:E, r e ) is a relativistic correction 

(|Challinor fc Lasenbvl 11998ft . which we take into 
account to first order in k^,T e /m e c 2 . The relativistic 
correction becomes significant when the electron tem- 
perature exceeds 10 keV, which is the regime of our 
cluster sample. 

One can combine Equations (HE]) and integrate along 
the line of sight to obtain the SZE in the apparent radi- 
ation intensity as 



A/. 



SZE 



j (i + e 2 /el 



,(l-3/3)/2 



(4) 



where 9 and 9 C are the angular equivalents of r and r c 
respectively. Because the clusters in our sample are not 
well resolved by AMiBA, we cannot get a good estimate 
of Iq, /?, and 9 C simultaneously from our data alone. In- 
stead, we use the X-ray derived values for /3 and r c , as 
summarized in Tabl e [H and then estimate the central 
specific intensity I (|Liu et al.ll2010ft by fitting Equation 
to the calibrated visibilities obtained by IWu et all 
). In the analysis we take into account the contam- 
ination from point sources and structures in the primary 
CMB. 

Given the /3-model described above, we can derive re- 
lations between cluster parameters and estimate them 
using the MCMC method. The parameters to be esti- 
mated are the electron temperature T e , r 2 5oo, total mass 
M t = Mt(r 2 5oo), gas mass M g = M g (r 250 o), and the 
integrated Compton Y = Yfasoa)- 

Theoretically M t (r 25 oo) can be fomula ted through 
the hydrostatic equilibrium e quation (e.g., iGrego et al.l 



tne n y 

unit] 



BonamcntEet_aL|[2008) 



M t (r 2 50o) 



3/?fc B T c 



' 2500 



G/im v 



(5) 



' 2500 



where G is the gravitational constant and /i is the mean 
mass per particle of gas in units of the mass of proton, 
m p . To calculate /i, we assume that /i takes the value 
ap propriate for cluster s with solar metallicity as given 
by lAnders & Grevessd (|1989f ). Here we use the value 
[i = 0.61. By combining Equation ([5]) and the definition 
°f r 2500i we can obt ain r 2 5oo as a function of /3, T , r c , 
and redshift z (e.g., iBonamente et a l. 2008) 



?*2500 



/3/3fc B T e 



(6) 



Giim p |7T/9 c (z) • 2500 

Then M g (r25oo) can be expressed, by integrating the 
c (r) in Equation ([I} as 



r/Dj 



2\ -3/9/2 



1 + T2 



l d9, 



M g (r) = ■i-n-fi e n e0 m p D 

.hi \ ", / 

where /i e = 1.17 is the mean particle mass per electron in 
unit of m p , Da is the angular diameter determined by z, 
and n e o is the central electron densit y, derived through 
the equation in LaRoque et al. (2006]): 



n e0 



AT m e c 2 r(|/3) 



f(x, T e )T C MBO-Tk B T e D A Tr^mf3 - \)9 C ' 



(8) 



where T is the gamma function, ATo is the SZE tempera- 
ture change, and Tcmb is the present CMB temperature. 
AT is derived as AT /T CM b = (e x - l)I /xe x I C mb. 

Finally, with the Iq computed earlier and the r 2 soo 
estimated here we can integrate the Compton y out to 
7"2500 to yield Y 



Y = 



2nAT Q 



f{x,T e )T C MB Jo 



9 2\ (l-3/3)/2 



9d9, (9) 



where 6^500 = ^500/^0 indicates the projected angular 
size of r 2500 . 

With the formulae as described above, for a set of /3, 
r c , and z as measured from X-ray observations and Iq 
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TABLE 1 

Parameters for isothermal spherical /3-model 









Without 100 kpc 


cut b 


With 100 kpc 


cut 6 


Cluster 


z 


D A 


/3 


r c 


A/ C 


P 




A/ C 






(Mpc) 




(") 


(xlO 5 Jy/sr) 




(") 


(xlO 5 Jy/sr) 


A1689 


0.183 


621 


609+ 005 
u.ou»_ 005 


26.6+° 7 7 


-3.13 ±0.95 


0.686t°;°l° 


48.0+j'^ 


-2.36 ±0.71 


A1995 


0.322 


948 


77n+ - 117 

u -' ' u -0.063 


38.9±ff 


-3.30 ± 1.17 


0-923±S:S3 




-3.19 ± 1.23 


A2142 


0.091 


340 


740+0-010 


188A±\l 2 2 


-2.09 ± 0.36 








A2163 


0.202 


672 


674+ - 011 


87.5±f| 


-3.24 ±0.56 


0.700+°,^ d 


78.8+^ d 


-3.64 ±0.61 


A2261 


0.224 


728 


°- 516 -0.013 


15.711;? 


-1.90 ±0.98 


U.b^S_ 020 


29.2+fl 


-2.59 ±0.90 


A2390 


0.232 


748 


600 +0 060a 


28.0l2;g a 


-2.04 ±0.65 


n crs+0.058 a 
u - oo -0.058 




-2.85 ±0.77 



a a 10% error is assumed for f3 and r c for which the original reference does not give an error estimation. 

b Referenc e -IReese et"aT1 j20u2fl for A1689, A1995, A2163, and A2261. ISanderson fc Ponmanl 020031 ); [Lancaster et~aT1 d2005ft 

for A2142. lAllen I J2000T) for A2390. 

c Best-fit values for AZq with foreground estimation from point sources and CMB IILiu et al.| [2010). 
d (3 fixed to a fiducial value 0.7 in[Boiiamcntc ct al. (2006), a 10% error is assumed. 

e Reference - IBonamente et al.l 1 120061) for A1689, A1995, A2163, and A2261. lAllen. et aT] (I2001D for A2390. 



from AMiBA SZE observation, we can arbitrarily assign 
a 'pseudo' electron temperature T c ^ y and then deter- 
mine the pseudo r25oo(T e (j)), M t (T oW ), M s (T e ^), and 
Y{T c{i) ). Given M t {T c{t) ) and M g (T e(i) ), we obtained the 
pseudo gas fraction f gas (T c(i )) = M g (T c ^)/M t {T c ^)). 
Using f g as{T e (i)) as a function of T e ^ we applied the 
MCMC method by varying T e and AIq to estimate the 
likelihood distribution of each cluster property. While es- 
timating the MCMC likelihood we assume that the like- 
lihoods of AIq and f gas are independent. The likelihood 
distributions of AIq for each clus ter are taken from the 
fitting results of iLiu et al.l (|2010D . while the likelihood 
distribution of f gas is assumed to be Gaussian with mean 
0.116 and standard deviation 0.005, which is the ensem- 
ble average ov er 38 clusters observe d by Chandra and 
OVRO/BIMA (jLaRoaue et al.ll2006h . 

In the process, the values of /?, r c , and z are taken 
from other observ ational results which are summarized in 
iKoch et all (|2010T ) and TableQ] We took the /3 model pa- 
rameters from both ROSAT and Chandra X-ray results. 
The Chandra results were derived by fitting an isother- 
mal /3 model to the X-ray data with a central 100-kpc 
cut. The aim of the cut-off is to exclude the complicated 
non-gravitational physics (e.g, radiative cooling and feed- 
back mechanisms) in cluster cores. Table [2] summarizes 
our results derived assuming an isothermal f3 model. We 
present the results obtained with isothermal f3 model 
parameters derived with and without 100-kpc cut both 
here. Figure [1] compares our results with the SZE-X- 
ray joint r esults obtained from OVRO/BIMA and Chan - 
dra data (jBonamente et al.l 120081 iMorandi et"alll2007f ). 
These are in good agreement. 

2.3. UTP P model 

The simulation done by [Hallman et al.l (|2007t) sug- 
gested incompatibility between isothermal /3 model pa- 
rameters fitted to X-ray surface brightness profiles and 
those fitted to SZE profiles. This incompatibility also 
causes bias in the estimates of Y and M g . They suggested 
a non-isothermal f3 model with a universal temperature 
profile (UTP). We also considered how the UTP j3 model 
changes our estimates of cluster properties in this section. 

In the UTP /3-model, the baryon density profile is the 



same as Equation fl]), and the t emperature profile can 
be written as (jHallman et al.ll2007t ): 

T e WHT) 500 T (l + (-^-) 2 ) , (10) 

where (T) 500 indicates the average spectral temperature 
inside r§QQ. Tq, a, and S are dimensionless parameters in 
the universal temperature profile model. S is the outer 
slope of the temperature profile, outside of a core with 
electron temperature T e0 = (T} 500 Tq. This core is of size 
otrsoQ. The total mass can be o btained by solving th e 
hydrostatic equilibrium equation (jFabricant et al.l ll980): 

In the isothermal /3-model, Equation (fTTj) can be reduced 
into the form of Equation However, in the UTP /3- 
model, the derivative of T e (r) with respect to r in Equa- 
tion (TTTT) is no longer zero. By applying Equation ([!} 
and Equation (fTU|) in Equation (jlll) . one can obtain: 

m, (r) . M± (4?rL + U + _£J) 

Gnm v \r z +r z c r 2 + a 2 r£ Q0 J\ a 2 r£ ao J 

By combining Equation (|12[) and the definition of r5oo, 
an analytical solution for rsoo can be obtained as: 

l(l + a 2 )(3l3A~r 2 c ) + 2SA + VD 

7-500 = v w^) ' (13) 

where A = 3k B T e0 (l + a' 2 )~ s /(4G/j,m p 7Tp c (z)-500), and 
D= [(l+a 2 )(^A-r 2 ) + 2SA} 2 + 8(l+a 2 )5Ar 2 . U8->0 
or a — > oo, which indicate the nearly isothermal case, 
Equation (fT3|) reduces to a form similar to Equation (|6|) . 

Using the definition of rsoo, M t (r 500) can be written 
as: 

M t (r 500 ) = 500 • -Ttrl ooPc (z). (14) 

For an arbitrary overdensity A, we can not find an an- 
alytical solution for arbitrary r& (i.e.: r25oo, ^200, etc.). 
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TABLE 2 

SZE DERIVED CLUSTER PROPERTIES IN ISOTHERMAL /? MODEL 



Cluster 



Without 100-kpc cut 
r-2500 k B T c M s M t 

(") (keV) (1O 13 M ) (1O 14 M ) 



Y 

(io- 10 ) 



»~2500 
(") 



k B T c 
(keV) 



With 100-kpc cut 

M g M t 
(1O 13 M ) (1O 14 M ) 



y 

(io- 10 ) 



A1689 
A1995 
A2142 
A2163 
A2261 
A2390 




214; 

159; 

237; 
172; 
174' 



-16 
-19 
-13 
-18 



-13 
-13 
-16 
-16 
-13 
-15 



10.0 
11 

15.4 
10.0 
11.9 



1.5 
1.6 
+ 1.7 
-2.3 



+ 1.6 
1.5 
+ 1.8 
1.7 
1.8 
1.9 



5.2 

8.5 

9.5 
4.6 
5.2 



1.3 

1.3 
+2.4 

2.5 



+1.6 
1.5 
1.3 
1.2 
+1.3 
1.2 



4.5 
7.5 

8.1; 
4.0; 

i.r 



1.2 
1.1 
+2.0 
2.3 



-1.5 
-1.3 
-1.2 
-1.0 
-1.1 
-1.1 



3.1 
1.9 



1.3 
1.2 
+1.0 
0.8 



-,+2.1 
-1.9 

2 2+ 11 
z ' z -0.9 
hi. 3 
-1.2 



8.0^ 



3.r 



10' 



10 



I I- far 

'' N= 3E 2 

| X 1 x =10.6229 

- 1 - dof=6 



10' 

kT (keV) 



HE 1 z =6.5637 

dof=4 



M (10 13 M ) 

g v sun' 



10 



10' 



T y 1 3 



Z =7.7302 
dof=5 



10' 



M,(10 14 M 



10' 



o 



10" 



10" 



Y(10" 



X =2.8344 
dof=5 



10' 



Fig. 1. — Comparison of T e (upper-left), M g (upper-right), Mt (lower-left), and Y (lower-right) of clusters derived from AMiBA 
SZE data based on iso t herma l /3 model with 100-kpc cut (x-axis) and t hose given in li t eratu re (y-axis). All y-axis values are 
from [Bonamcntc ct al. (2008), exc ept for the Y values , which are from Morandi ct al. (200 ^), and those for A 2390, which is 
indicated by a circle with T e from IBenson et all (|2004D and Mt calculated from the data in lBenson et al.l l|2004l ). The dashed 
lines indicate y = x. 



However, with the known r§Qo, we can still find the nu- 
merical solution for r& easily. We can then solve for 
M t {r&) using Equation ([12]) , 

To yield the central electron number density, we con- 
sider the formula for the Compton y resulting from 
the U TP /3-model (see the Appendix of Hal lman et al.l 
( 2007)). By s etting the proj e cted r adius b = in Equa- 
tion (A10) in Hallman et al.l ()2007f ). one can obtain: 



n e0 



AT m e c 2 



f(x, T e )T C MBcr T k B (T} 500 T o I sz (0) 



(15) 



where 



isz(0) = 



r(f + 6) 

and F2.1 is Gauss' hypergeometric function. Here we 
assume f(x,T e ) = f(x, (T} 500 Tq), and the change of 
f(x,T e ) due to the change of T e along line of sight is 
negligible. Actually, by numerical calculation we found 
that the error in Equation (1151) caused by this assump- 
tion is less than 1%. Because the UTP j3 model assumes 
the electron density profile as same as the isothermal (3 
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model, we can rewrite M g in UTP model by simply ap- 
plying Equation (fT5|) in Equation 

Thus, the integration of the Compton y profile, instead 
of Equation ((9|), becomes: 



Y = Y Q 



1 + T2 



(l-3/3)/2 



y 500 



-6 



I ,/.,, ! F{6)9d9, 
(17) 

where Y = (2^ATo)/(/T C mb-F(0)) and F(0) = 
F 2A (S, 1/2; 3/3/2 + 6,1- (r 2 + 9 2 )/(a 2 r 2 00 + 2 )). 

We were not able to constrain the parameters f3, r c , 
S, and a of the UTP significantl y with our SZE data 
alone. However, the simulation of lHallman "eTall (|2007h 
suggested that there is no significant systematic differ- 
ence between the values of /3 and r c resulting from fit- 
ting an isothermal (3 model to mock X-ray observations 
and those parameters fitted using the UTP (3 model. 
Therefore, we simply assume that the ratio between the 
isothermal /3j SO value and UTP /3utp value is 1 ± 0.1, 
and r C! j SO /r Ci uTP = 1 ± 0.2, for each cluster. We also 
assume 5 — 0.5, a = 1, and To = 1.3. Those values 
are ta ken from the average of results of Hallman et al. 
(2007). Then we fit A/ to AMiBA SZE observation 
data with the UTP j3 model parameters above by fix- 
ing 5, a, and To, and treating the likelihood distribu- 
tions of /?utp and r Cj uTP as two independent Gaussian- 
distributions. Finally, we applied the MCMC method, 
which varies Aio, r c , and (T) 500 , to estimate cluster 
properties with the equations derived from the UTP f3 
model and the data fitting results. 

Table [3] summarizes our results derived with the UTP 
j3 model. Figure [5] compares our results with the SZE-X- 
ray joint r esults obtained from OVRO/BIMA and Chan - 
dra data (jBonamente et all 120081 : IMorandi et"aD l2007f l. 
These are also in good agreement. We find that the 
electron temperature derived with the UTP j3 model are 
in significantly better agreement with the temperatures 
from Chandra X-ray measurements. 

3. EMBEDDED SCALING RELATIONS 



The self-similar model (Kaiser] ri986| ) predicts simple 
power -law scaling relations between cluster prop erties 
(e.g., iBonamente et al.ll2008l: IMorandi et al.ll2007l ). Mo- 
tivated by this, people usually investigate the scaling re- 
lations between the derived cluster properties from ob- 
servational data to see whether they are consistent with 
the self-similar model. However, the method described 
above is based on the isothermal /3-model and the UTP 
/3-model. Therefore, there could be some embedded re- 
lations which agree with self-similar model predictions 
between the derived properties. We investigated the em- 
bedded relations through both analytical and numerical 
methods. 

3.1. Analytical formalism and numerical analysis 

In the isothermal j3 model, by applying Equation ([5]) 
in Equation ([5]), M t can be rewritten as 

M t = 2500 X Pc (z) ("P^— -J— n 
3 \ (j\im v 2500 • |7r/9 c (z) 

(18) 

A_s we can see, while f3 is set to be a constant, and 
which implies 3(3kBT /(Gfj,m p ■ 2500 • 



f 7rp c (z)) >> r 2 , the relation M t ocTe will be obtained. 
However, for some of the clusters we considered in this 
paper, the values of r25oo/ r c are only slightly above 2. 
Therefore, we have to investigate the scaling relation be- 
tween M t and T c by considering <91nM t /91nT . 

By partially differentiating Equation (|18|) by T c , and 
multiplying it by T e /M%, we can get 



d In Af t = 3 (r 2 500 + r c 2 
91nT e 2 



(19) 



' 2500 



which decreases from 1.875 at r25ooAc = 2 to 1.5 as 
^2500 Ac — > co. That implies M t behaves as M t <x T e 1875 
while r 2 5ooAc ~ 2 and M t oc T^ 5 while r 2 5ooAc ap- 
proaches infinity. This result shows that there is an 
embedded M t -T e relation consistent with the self-similar 
model in the method described above. 

If we assume that the gas fraction / gas is a constant, 
the scaling relation between M g and T c will be as same 
as the relation between M t and T c . 

In order to investigate the relations between integrated 
Y and the other cluster properties, we consider Equa- 
tion ©. By combining Equation (JS])-©, one can obtain: 



AT = 



M g (r 2500 )f{x, T e )T CMB G T k B T c T (§0 - \) 



^/^ c m p Dim cC *T (|/3) J^ 2500 ( 1 + 



5? 



-3/3/2 



9 2 d6 



(20) 

Then we combine Equation (|2"0"|) and Equation ^ and 
obtain: 



Y = 



Tr 1/2 M g {r 250Q )<j T k B T c 
2^ e m v m e c 2 D\ 



'(^2500)^c>^)) 



(21) 



where 



ff(#2500, 6c, P) — 



r m j:- (i + g) 



(l-3/3)/2 



-3/3/2 



2 d6 



' 2500 



>> rf 



( 22 ) 

is a dimensionless function of ^2500: and j3. 

We also calculated d In Y/d In T c to investigate the be- 
havior of Y when T e varies (see Figure [3]). As we can see 
in Figure [31 d\nY/d\nT c varies between 2.45 and 2.75 
while r 2 5ooAc > 2.0 and 0.5 < /3 < 1.2. We also noticed 
that 9 In Y"/ <9 In approaches 2.5 as 72500 Ac approaches 
infinity. This result indicates that behaviour similar to 
the self-similar model is built into scaling relation studies 
based solely on SZE data. 

The effect of varying j3 is investigated. If we consider 
power law scaling relation 

Q = W A X B (23) 

between M t and T e with M t written as Equation (fTS)) . 
one can find that changing the value of j3 will only affect 
the normalization factor A. In other words, if we change 
/3 to 13', A will be changed to A' = A + B log 10 (/3'//3). 

In the Y-T e relation, f3 will affect the scaling power B 
as shown in Figure [3] B varies within a range of only 
0.04 while 0.5 < (3 < 1.2. 

Considering the UTP (3 model, we undertook a similar 
analysis of the embedded scaling relation. The results, 
which are similar with those obtained with the isothermal 
j3 model, are shown in Figure HJ 
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TABLE 3 

SZE DERIVED CLUSTER PROPERTIES IN THE UTP MODEL 







Without 100-kpc cut 








With 100-kpc cut 




Cluster 


^2500 




M g 


M t 


y 


^2500 
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(io- 10 ) 
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458l« 
245l?| 
160ig 
16611? 


8.9ti : g 

9.711;? 
13.0±1;S 
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8.i±l:J 


°- D -1.8 

7.6±1;| 
7.6+11 
10.1+i? 
3.6±i| 

4.5+1:1 


4.8li;I 

7+2.3 
D " ' -2.3 

6.41 2 ;! 
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».8_ 2 4 

3.9±1;I 


4+1.7 

1 8+ 1 ' 1 
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1.518:1 
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a The average electron temperature up to rsoo (i- e - : {^)5oo ^ n Equation dl0| l) 
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Fig. 2. — Comparison of T e (upper-left), M g (upper-right), M t (lower-left), and Y (lower-right) of clusters derived from AMiBA 
SZE d ata based on the UTP j3 model with 100-kpc cut (x-axis) and t hose given in literatur e (y-axis). All y-axis values are 
from iBonamcn te et aD (|2008l ) , exc ept for the Y values , which are from Morandi et al.l (|2007l), and those for A 2390, which is 
indicated by a circle with T e from Be nson et al.l (|2004D and M t calculated from the data in Benson et al. (20041). The dashed 
lines indicate y = x. 



3.2. Calculation oj Scaling Relations 

Here we investigate the Y-T c , Y-M t , and Y-M g scaling 
relations for the quantities derived above. We also study 
the M t -T c scaling relation with t he M t from AMiBA SZE 
data and the r e from X-ray data (Bonam ente et ai]|2008l : 
iMorandi et aLl[2007l ). 

For a pair of cluster properties Q-X, we consider the 
power-law scaling relation (Equation ([2"3)) ). To estimate 
A and B, we perform a maximum-likelihood analysis 
in the log-log plane. For the M t -T e relation, because 
M t and T e are independent measurements from differ- 



ent observational data, we can simply per form linear 
minimum- y 2 analysis to estimate A and B (| Press et al.l 
119921: iBenson et atl 120041) . On the other hand, for the 
SZE-derived properties, because they are correlated and 
so are their likelihoods (i.e., L(Q,X) ^ L(Q)L(X), as 
manifested by the colored areas in Figure 0, we cannot 
apply x 2 analysis. Instead we use a Monte Carlo method 
by randomly choosing one MCMC iteration from each 
cluster many times. With each set of iterations we de- 
rived a pair of Ai and Bi using linear regression method. 
Finally we estimate the likelihood distribution of A and 
B using the distribution of {Ai} and {Bi}. The results 
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Wc 

Fig. 3. — Embedded scaling relation between Y and T e . The 
shaded scale indicates different /3 from 0.5 (the darkest line) to 
1.2 (the lightest line). The dashed line indicates the predicted 
value by self-similar model. 
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Fig. 4. — Embedded M t -T c (upper panel) and Y-T c (lower 
panel) scaling relations in UTP /3 model. The grey scales 
indicate different ft from 0.5 (the darkest line) to 1.2 (the 
lightest line). The dashed lines indicate the predicted values 
by self-similar model. 

are presented in Table 0] and Figures [5] and [6j However, 
as we discussed in Section 13- H the scaling relations be- 
tween SZE-derived properties should be interpreted as 
a test of embedded scaling relations rather than estima- 
tions of the true scaling relations. On the other hand, 
the M t -T e relation compared M t and T e from different 
experiments. Therefore, we can regard it as a test of the 
scaling relation prediction. 

4. DISCUSSIONS AND CONCLUSION 

We derived the cluster properties, including T c , 
r 2500 ,M t , Mg and Y, for six massive galaxy clusters 
(M t (r 25 oo) > 2 x 10 14 M Q ) mainly based on the AMiBA 
SZE data. These results are in good agreement with 
those obtained solely from the OVRO/BIMA SZE data, 
and those from the joint SZE-X-ray analysis of Chandra- 
OVRO/BIMA data. In the comparison, the SZE-X-ray 
joint analysis gives smaller error bars than the pure SZE 
results, because currently the uncertainty in the mea- 
surement of the SZE flux is still large. On the other 
hand, in our current SZE-based analysis, due to the in- 
sufficient u-v coverage of the 7-element AMiBA we still 



TABLE 4 

Scaling relations of SZE-derived cluster 
properties 



Scaling 
Relations 


A 


B 


B thy 


D\E(z)Y,T 
D 2 A E(z)- 2 / 3 Y,M t 
D 2 A E(z)- 2 / 3 Y, M g 
E(z)M t ,T 


-4.321°;^ 

-4.80l°; 2 21 
-4 89+ - 22 

<±.o»_ 22 
°- 66 -0.12 


2.48+8:3° 

1 29+ 028 

n 05+O.66 
U.yO_o.60 


2.50 

1.67 
1.67 
1.50 



Note. — All cluster properties used in the analysis 
are based on the AMiBA SZE data (see Sec. [2j, ex- 
cept for the T in the M-T relation, where the T is 
from IBonamente ct al. I 2008 f) for A168 9, A1995, A2163, 
A2261, and from Morandi c t al.l ( 1200711 for A2390. The 
units of T, D 2 A Y, M t , and M g are 7keV, Mpc 2 , 1O 14 M , 
and 1O 13 M0, respectively. The last column Bthy in- 
dicates the theoretical values predicted by self-similar 
model. In the first column. E 2 (z) ~ £1m (1 + z ) 3 + 

(i-n M -n A )(i + z) 2 + n A .' 




0.2 0.4 0.6 0.8 1 1.2 1.4 

log(M l /10 14 M sun ) 




Fig. 5. — Scaling Relations of Y-T e (upper), Y-M g (middle), 
and Y-Mt (lower) based on the AMiBA SZE derived results. 
Gray areas indicate the 68% confidence regions for the pa- 
rameter pairs of each cluster. Solid lines are the best fits as 
in Tab. g| 

need to use X-ray parameters for the cluster model i.e., 
the (3 and 9 C for the /3-model. However, iNord et"aLl 
(2009) have deduced /3 and 9 C from an APEX SZE ob- 
servation alone recently. For AMiBA, the situation will 
be improved when it expands to its 13- element configu - 
ration with 1.2m antennas (AMiBA13; IHo et al.|[2009() . 
and thus much stronger constraints on the cluster proper- 
ties than current AMiBA results are expected. Further- 
more, with about three times higher angular solution, 
we should be able to estimate /3 and 9 C from our SZE 
data with AM1BA13 and make our analy sis purely SZE 
based (|Ho et al.1120091: iMolnar et al.ll2010[ ). Nevertheless, 
the techniques of using SZE data solely to estimate clus- 
ter properties are still important, because many upcom- 
ing SZE surveys will observe SZE clusters for which no 
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Fig. 6. — M t — T e scaling; relation between the X-ray mea- 
sured T e (|Bonamente et all I2008I ; IMorandi et~aTl 120071 ) and 
the AMiBA derived M t . The boxes indicate the la errors for 
each cluster. The solid line is the best fit as in Tab. [3] 

X-ray data a re available (iRuhl et al.l l200l lFowlerll2004 
iKanek o 2006; Ho et al.ll2009f ). especially for those at high 
r edshifts. 

lHallman et al.l |2007) suggested that adopting the 
UTP /3 model for SZE data on galaxy clusters will reduce 
the overestimation of the integrated Compton Y500 and 
gas mass. However, the Y2500 values we obtained with the 
UTP model are not smaller than those obtained with the 
isothermal model. The M g (r25oo) values deduced using 
the UTP model are even larger than those deduced using 
the isothermal model. 

For the case of integrated Compton Y, when we com- 
pare Y500 deduced using the UTP model Ysoo.utp, and 
those deduced using the isothermal model Ysoo.iso, we 
found tha t the Ysqo.utp are sm aller than l5oo,iso, as pre- 
dicted by lHallman et al.l (|2007f ) . The reason is that the 
Compton y profile predicted using the UTP (3 model will 
decrease more quickly than the profile predicted by the 
isothermal f3 model, with increasing radius. Therefore, 
the ratio Ya^utp/Yaaso will decrease as A decreases. 

We also noticed that the electron temperature val- 
ues obtained with the isothermal model are significantly 
higher than the temperatures deduced from X-ray data 
for most clusters we considered. The temperatures of 
clusters obtained using the UTP model are lower than 
those obtained with the isothermal model and thus are 
in better agreement with those deduced from X-ray data. 
Therefore, in the UTP model, with similar Y2500 and 
lower temperature, we should get larger M g . 

The electron temperatures derived using the UTP j3 
model are in better agreement with X-ray observation 
results than those derived using the isothermal (3 model. 
This result implies that the UTP /3 model may provide 
better estimates of the electron temperature when we 
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can use only the f3 model parameters from X-ray obser- 
vation. However, we noticed that the UTP f3 model pro- 
duced larger errorbars than the isothermal f3 model did. 
These increased errors are based on the uncertainties of 
(3 and r c which we insert by hand. On the other hand, 
because we treat /3 and r c as independent parameters in 
this work, the uncertainty could be over estimated due to 
the degeneracy between these two parameters. If we can 
access to the likelihood distributions of f3 and r c of the 
UTP f3 model derived from observation, the error-bars 
might be reduced significantly. 

There is a concern that the scaling relations among the 
purely SZE-derived cluster properties may be implicitly 
embedded in the formalism we used here. In this pa- 
per, we also investigate for the first time the embedded 
scaling relations between the SZE-derived cluster prop- 
erties. Our analytical and numerical analyses both sug- 
gest that there are embedded scaling relations between 
SZE-derived cluster properties, with both the isothermal 
model and the UTP model, while we fix /?. The embed- 
ded Y-T and M-T scaling relations are close to the pre- 
dictions of self-similar model. The results imply that the 
assumptions built in the pure-SZE method significantly 
affect the scaling relation between the SZE-derived prop- 
erties. Therefore, we should treat those scaling relations 
carefully 

Our results suggest the possibility of measuring cluster 
parameters with SZE observation alone. The agreement 
between our results and those from the literature pro- 
vides not only confidence for our project but also sup- 
ports our understanding of galaxy clusters. The upcom- 
ing expanded AMiBA with higher sensitivity and better 
resolution will significantly improve the constraints on 
these cluster properties. In addition, an improved de- 
termination of the u-v space structure of the clusters 
directly from AMiBA will make it possible to measure 
the properties of clusters which currently do not have 
good X-ray data. The ability to estimate cluster proper- 
ties based on SZE data will improve the study of mass 
distribution at high redshifts. On the other hand, the 
fact that the assumptions of cluster mass and tempera- 
ture profiles significantly bias the estimations of scaling 
relations should be also noticed and treated carefully 
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